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and  ultimately  cost  of  many  circuits.  We  wish  to  demonstrate  the  utility  of  Raman  spectroscopy  as  a  tool 
for  the  development  and  design  of  silicon  microstructures. 

The  equations  for  a  general  two-dimensional  stress  field  are  discussed.  Calibration  studies  using 
macromechanical  fixtures  for  single  crystal  silicon  specimens  under  a  two-dimensional  stress  field  are 
presented.  Our  measurements  show  good  agreement  with  the  theoretical  values  and  thus  validate  the 
approach  taken.  Stress  maps  of  conventionally  fabricated  test  structures,  laser  machined  structures,  and 
polysilicon  structures  are  presented. 
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1.  INTRODUCTION 


High  stress  levels  play  a  significant  role  in  determining  the  life,  the  reliability,  and  the  failures  of 
micromachined  and  microelectronic  devices.  High  stresses  due  to  coefficient  of  thermal  expansion 
(CTE)  mismatch  and  discontinuities  have  led  to  failures  in  interconnects  and  metallizations  due  to  stress 
migration  and  electromigration  effects.  Stress  in  films  may  also  cause  adhesion  failure  between  film 
layers.  Fatigue  life  may  be  considerably  shortened  by  the  presence  of  localized  high  stress  in  critical 
regions.  Structural  designs,  workmanship,  and  manufacturing  methods  including  laser  machining  also 
play  a  role  in  determining  stress.  We  have  proposed  Raman  spectroscopy  as  a  noncontact  method  for 
measuring  highly  localized  stresses  at  the  surface  of  silicon  structures  to  understand  the  effects  of  macro- 
and  micromaterial  property  variations,  design,  and  defects  as  they  influence  stress  concentrations. 

Raman  spectroscopy  has  excellent  spatial  resolution  of  1  |am  and  is  well  suited  to  probing  the 
nanostructures  used  in  MEMS  devices.  The  stress  measurement  is  based  on  the  relative  frequency  shift  of 
the  Raman  spectra  when  the  crystal  lattice  is  strained.  Stresses  as  small  as  8.7  MPa  (1.25  ksi)  may  be 
detected. 

This  paper  presents  the  Raman  secular  equation  for  cubic  crystal  systems  and  discusses  the  general,  two- 
dimensional  surface  stress  field  effects  on  Raman  fi'equency  shifts  and  polarization.  Raman  stress 
measurements  of  single-crystal  silicon  test  structures  undergoing  four-point  bend  loading  are  used  to 
validate  the  theory.  These  are  the  first  measurements  reported  for  single  crystal  silicon  undergoing 
tensile  stress.  Residual  stress  measurements  in  polysilicon  cantilever  beams  and  preliminary  stress  maps 
in  laser  machined  test  structures  using  the  Revise  Inc.  (see  Reference  1)  process  are  also  reported. 
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2.  THE  RAMAN  EQUATIONS  FOR  CUBIC  SYSTEMS 


The  Raman  dynamic  equations  for  optical  modes  in  a  cubic  crystal  were  used  to  derive  the  stress 
dependent  frequency  shift  relationships.  The  secular  equation  for  a  cubic  crystal  such  as  silicon  is 
expressed 


[pe  -fq(e  +e  )-A]  2re 
•■^xx  yy  zz^  xy 


xz 


2ie 


xy 


[pe  +q(e  +e  )-A] 
yy  XX 


2rc 


2ie 


xz 


2re 


yz 


[pe  +q(e  +e  )'-A] 
zz  XX  yy^ 


=  0 


(1) 


where  p,  q,  and  r  relate  to  the  atomic  spring  constants  for  the  crystal,  the  en  are  the  strain  components  and 
X  =  Q  -  COo  represents  the  frequency  shift.  Since  X  is  so  much  smaller  than  cOo,  the  shifted  frequency  Q. 
can  be  approximated  by  £2  =  cOq  +  X/2(Oo  and  A£2  =  O  -  ©o  =  X/2(Oo. 

The  two-dimensional  stress  field  is  appropriate  to  a  surface  where  stress  perpendicular  to  the  surface  is 
absent.  Consider  a  general  two-dimensional  stress  field  Xn,  T22,  and  X12  where  Xu  and  X22  are  two  normal 
stresses  in  directions  1  and  2,  respectively,  and  X12  is  the  shear  stress  with  respect  to  a  Cartesian 
coordinate  system  with  axes  1,  2,  and  3.  Further,  let  the  1,  2,  and  3  axes  coincide  with  the  principal  axes 
of  the  cubic  crystal  material  with  axis  3  normal  to  the  surface.  Let  the  three  stress  components  be  related 
by  the  parameters  a  and  P  as  follows: 

'Til  =^1,  1^22  =  ax,  and  X12  =  P  x  (2) 

The  following  three  roots  (frequency  shifts)  are  derived  from  solving  Eqn.  (1): 

A£2,  =  (x/4©o){p(l+a)(S„  +  S12)  +  q(l+a)  (S,,  +  3S,2)  +  [(l-a)'(Su  -  Si2)'(p-q)' 

+4(rS44P)']'''}, 

AQ2=  (x/4Wo){p(l+a)(S„  +  S,2)  +  q(l+a)  (Sn  -t  3Si2)  -  [(l-a)"(S„  -  S,2)"(p-q)" 

+4(rS44p)']"''}, 

AQ3  =  (x/2©o){p(l+a)Si2  +  q(l+a)  (Sn  +  S12)}  (3) 

where  Sy  are  elements  of  the  compliance  matrix  that  relates  stress  to  strain;  p,  q,  and  r  are  calibration 
constants;  and  ©0  is  the  reference  frequency  of  unstressed  silicon,  520  cm"'.  It  is  interesting  to  note  that 
the  root  AQ3  is  independent  of  the  shear  stress  whereas  the  other  two  roots,  AQi  and  AQ2  will  show  a 
splitting  relative  to  each  other  that  is  shear  and  biaxial  stress  dependent.  The  values  of  p/  ©0^  =  -1.85, 
q/  ©o^  =  -2.31,  r/  ©0^  =  -0.71,  S„=  7.68  x  10’'^  Pa‘‘,  812=  -2.14  x  10''^  Pa'*,  and  844=  12.7  x  10'*^  Pa'*  were 
used  for  single  crystal  silicon  (see  Reference  2).  For  polycrystalline  silicon  (see  Reference  3)  the  values 
are  p/  ©0^  =  -1.40,  q/  ©0^  =  -2.00,  r/  ©0^  =  -0.67,  8,1=  7.73  x  lO'*^  Pa'*,  812=  -2.15  x  10'*^  Pa'*,  and  844= 

12.7  X  10'*^  Pa'*. 


2.1  Transition  Dipole  Moments  and  Selection  Rules 

The  three  frequency  shifts  associated  with  the  stress  field  are  not  necessarily  observable  simultaneously 
at  a  specific  orientation.  Their  detectability  depends  on  the  combination  of  the  polarization  of  the 
incident  light  and  the  scattered  light.  For  cubic  crystals  such  as  silicon,  the  Raman  tensors  in  the  1, 2,  and 
3  directions  are  respectively 
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The  Raman  tensors  corresponding  to  the  three  Eigenvectors  under  the  above  two-dimensional  stress  field 
are 

0  d  0 

d  0  0 

.0  0  0. 

where  d  is  the  lattice  constant.  Its  magnitude  is  of  the  order  of  10  nm.  Ai  and  Aa  are  defined  in  Eqn.  (5d) 

Selection  rules  were  examined  by  evaluating  the  dipole  moment  for  the  strain-induced  Eigenvectors  for 
the  following  polarizations  of  incident  and  scattered  light.  Let  Ci  and  Cs  be  the  incident  and  scattered  light 
with  polarizations  of: 


ll. 

Ik 

e]  = 
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- 1 

_ 1 

where  lim  are  directional  cosines  and  +  (1*2)^ +(^*3)^  *=1  k  =  i,s. 
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The  superscript]  labels  the  three  roots  in  (3).  The  scattering  intensity  of  the  Raman  signal  is  defined  as 
C2  |e*i  Rj  where  C  is  a  scalar  constant.  The  results  summarized  below  show  the  amplitude  strength 
of  the  observed  frequency  shift  for  each  solution  in  (3). 

e‘  R el  =  i;, (Z;, At  +ll, )d / yjl  +  A],  (5a) 

ef  Rcl  =  lf,il^,+l^,A2)d/yll  +  A\,  (5b) 

ef  Rel  =  (5c) 

where 

.  (l-a)(Sn-Siz)(p-q)/2-D  . _ _ 

rSujS  ’  (,l-a)(Sn-Si2)ip-q)/2  +  D' 

(5d) 

D=i4(~l  -a)‘^iSu -Si2fip-qf  +  4irSu/3f . 


A  strategy  to  uniquely  identify  the  transitions  from  <001>  single  crystal  silicon  can  be  developed  from 
the  selection  rules.  It  is  seen  from  Eqn.  (5)  that  for  cases  corresponding  to  backscatter  from  a  <001?» 
surface,  only  longitudinal  phonons  can  be  observed  since  only  (5c)  remains  nonzero  for  polarizations  in 
the  x-y  plane.  When  the  polarization  angles  of  incident  light  and  the  scattered  light  complement  each 
other  in  backscatter  from  the  <00 1>  surface,  the  amplitude  (5c)  is  maximized.  When  the  incident  and 
scattered  polarizations  are  parallel  (^1,  =  0)  and  the  polarization  vectors  are  aligned  with  either  the  x  or  y 
axis  in  the  backscatter  configuration,  all  amplitudes  become  zero.  The  (5c)  term  is  nonzero  for  all  other 
polarization  orientations.  Hence  its  behavior  can  be  used  to  determine  the  crystallographic  orientation  of 
the  cubic  crystal,  (see  Reference  4). 
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In  order  to  observe  the  effects  of  shear  as  expressed  by  AQ]  and  AQ2.  the  sample  must  be  tilted  to  allow  a 
nonzero  projection  along  the  z-axis  to  allow  non-vanishing  or  (.], .  In  this  tilted  configuration,  if  the 
polarizations  of  the  incident  and  scattered  light  are  parallel,  then  the  longitudinal  phonon  transition  (5c) 
can  be  nulled  to  enable  observation  of  only  the  first  two  roots.  In  general,  as  observed  in  the  backward 
direction  to  the  tilted  incident  beam,  one  will  detect  the  sum  of  the  transition  intensities  of  (5a)  and  (5b). 
If  the  frequency  difference  between  AQ(j=i)and  AQ(j=2)is  sufficiently  large  compared  to  the  linewidth, 
then  it  will  be  very  easy  to  identify  each  transition  by  a  rotation  of  the  polarization  vector  around  the 
incident  axis.  One  or  the  other  transition  intensity  will  go  to  zero  at  an  angle  determined  by  Aj.  In 
general,  the  ease  of  identifying  the  transitions  and  shifts  is  dependent  on  how  muc^iaxial-  and  shear- 
stress  is  present.  The  separation  of  the  peak  positions  of  AQ|  and  AQ2  is  equal  to  —  where  D  is  defined 

above.  In  principle,  if  there  is  knowledge  of  the  frequency  shifts  of  all  three  transitions  and  the 
polarization  angle  at  which  either  (5a)  or  (5b)  becomes  zero,  then  it  is  possible  to  derive  the  values  of  x, 
a,  and  p. 
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3.  EXPERIMENTAL  SETUP 


In  the  Raman  stress  experimental  setup,  a  monochromatic  light  or  laser  shines  on  the  surface  of  a 
specimen  through  a  microscope.  The  specimen  sits  on  a  translation  stage  whose  motion  is  controlled 
through  a  micrometer.  The  Raman  scattered  light  is  analyzed  with  a  spectrometer.  When  the  specimen  is 
unstressed,  the  spectrometer  measures  a  reference  spectrum  whose  center  position  is  that  of  the 
unstressed  sample.  When  the  specimen  is  placed  in  a  stressed  state,  the  spectrum  will  show  a  frequency 
shift  that  varies  linearly  with  the  applied  stress.  The  measurements  are  conducted  using  an  argon  ion 
laser  with  a  wavelength  of  488  nm.  The  diameter  of  the  laser  spot  is  typically  1  to  2  p,m  with  an  estimated 
penetration  depth  for  the  incident  light  of  0.6  pm.  The  Raman  spectra  were  least-square  fitted  with  a 
Lorentzian  profile  to  determine  the  center  frequency  shift. 

The  specimens  that  were  used  in  the  experimental  tests  included  both  calibration  and  MEMS  specimens. 
The  calibration  specimen  was  a  piece  of  single  crystal  silicon  in  a  strip  configuration.  The  strip  specimen 
was  either  clamped  or  simply  supported  at  its  two  end  regions.  The  specimens  had  a  clear  central  span  L 
and  were  loaded  with  two  rollers  located  symmetrically  with  respect  to  the  center  of  the  specimen.  Loads 
in  the  samples  were  induced  by  turning  the  set  screw  located  at  the  center  of  the  supported  strip.  The 
magnitude  of  the  load  setting  was  recorded  from  the  load  cell.  The  maximum  surface  stress  of  the  sample 
between  the  central  supports  is: 

CT,  =  3PbV(Lwt^)  (6) 

for  the  clamped  sample  and  for  the  simply  supported  sample  is: 

CTs  =  3Pb/(wt^)  (7) 

where  P  is  the  total  load  measured  by  the  load  cell,  b  and  L  are  the  dimensions  shown  in  Figure  1,  and  w 
and  t  are  the  width  and  thickness,  respectively,  of  the  strip. 

Strip  samples  of  <001>  single  crystal  silicon  were  fabricated  from  3-inch  diameter  standard  n-type  single 
crystal  silicon  wafers  for  use  in  either  fixture.  The  sample  dimensions  were  12.7-mm  (0.500  in.)  in  width 
and  0.394-mm  (0.0155  in.)  in  thickness.  They  consist  of  two  types  such  that  a  unidirectional  stress  could 
be  applied  either  in  the  (1 10)  or  (100)  directions. 

The  calculated  results  for  single  crystal  silicon  using  (2)  relate  the  Raman  shifts  in  units  of  cm’’  to  stress 
a  in  units  of  Pascals  as  follows: 

For  applied  uniaxial  stress  along  the  (100)  direction,  x  =  as=  a,  a  =  0  and  P  =  0:AQ3  =  -2.3  x  lO'^Oj. 

For  applied  uniaxial  stress  along  the  (1 10)  direction,  x  =  0^/2=  a/2,  a  =  1  and  P  =  1:  AQ3  =  -2.3 

X  lO  ’aj. 

The  calculated  results  for  polysilicon  using  (2)  relate  the  Raman  shifts  in  units  of  cm'*  to  stress  x  in  units 
of  Pascals  as  follows: 

For  “normal”  uniaxial  stress,  a  =  0  and  P  =  0:  A£23  =  -2. 12  x  10  ®x. 

For  “normal”  equal  biaxial  stress  without  shear,  a  =  1  and  P  =  0:  A£23  =  -4.24  x  lO’^x. 

The  term  “normal”  refers  to  stress  components  in  the  same  direction  as  the  crystal  axis.  The  latter  case 
may  correspond  to  stress  caused  by  a  mismatch  in  the  coefficient  of  thermal  expansion  between  two 
adjacent  layers. 

The  MEMS  specimen  is  a  cantilever  beam  fabricated  by  the  MUMPS  process  (see  reference  5).  Its 
nominal  dimensions  are  200p.m  long,  20pm  wide,  and  2pm  thick.  It  simulates  an  electrical  relay.  Since 
polysilicon  is  partially  conductive,  an  insulating  oxidation  layer  was  deposited  between  the  cantilever 
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and  the  ground  to  allow  application  of  a  voltage  that  induces  an  electrostatic  force  in  the  cantilever. 

When  the  voltage  increases  to  a  critical  value,  the  cantilever  beam  will  collapse  onto  the  ground.  The 
deformation  induced  by  the  voltage  will  generate  bending  stresses  in  the  beam  with  maximum  values  at 
the  clamped  end.  Figure  2  shows  a  photograph  of  a  set  of  three  cantilever  beams. 

Simple  test  structures  were  fabricated  from  <100>  single  crystal  silicon  using  a  Revise  Inc.  (see  reference 
6)  laser  machining  work  station.  The  process  uses  thermal  heating  by  a  laser  to  enhance  chemical  etch 
rates  in  an  atmosphere  of  chlorine  gas.  A  resolution  of  2  jim  is  possible  when  the  laser  beam  is  focused 
well.  Square-shaped  trenches  were  fabricated  by  raster  scanning  the  laser  spot  at  a  power  of  a  few  watts 
with  the  5 14.5  and  488.0-nm  wavelengths  from  an  Ar  ion  laser. 


Figure  1.  Schematics  of  a  strip  specimen  in  a  Figure  2.  Photograph  of  MEMS  cantilever  beams 

fixture  for  either  a  clamped  or  simply  fabricated  by  MUMPS  process, 

supported  end. 
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4.  RESULTS  AND  DISCUSSION 


Figure  3  depicts  the  Raman  spectra  of  wafer  cut  from  an  annealed  single  crystal  silicon  boule.  The 
average  center  frequency  shift  measured  at  random  locations  on  the  wafer  is  520.28  +  0.06-cm''  at  room 
temperature.  This  uncertainty  should  correspond  only  to  a  uniaxial  stress  with  no  shear  component  with  a 
value  of  +  26-MPa  (or  +3.8-ksi).  The  estimated  error  in  the  determination  of  the  fitted  line  position  is 
0.02-cm''  or  8.7-Mpa  (1.25-ksi) .  Measurements  were  also  taken  on  an  annealed  polysilicon  layer 
deposited  on  a  SiO:  layer  on  top  of  a  single-crystal  silicon  substrate  in  Figure  4.  These  measurements 
showed  a  spatial  variability  of  +  68.3-MPa  (9.8-ksi)  which  is  well  in  excess  of  the  residual  stress  of  the 
uncoated  substrate,  indicating  that  the  coating  has  increased  the  residual  stress.  This  finding  is  consistent 
with  the  earlier  work  of  Miura  (see  Reference  7). 

Raman  frequency  measurements  for  <001>  single  crystal  silicon  strip  with  (100)  applied  unidirectional 
stress  are  presented  in  Figure  5.  The  measured  slope  is  2.86  +  0.27  x  10'^-cm''/Pa.  The  theoretical  slope 
is  2.3  X  10'^'cm''/Pa.  The  Raman  frequency  measurements  with  (110)  applied  unidirectional  stress  are 
presented  in  Figure  6.  This  sample  and  its  stress  orientation  correspond  to  the  equal  biaxial  case 
discussed  above.  The  measured  slope  is  2.07  +  0.14  x  10'’-cm''/Pa  compared  to  the  expected  value  of  2.3 
X  10'®'cm’'/Pa.  The  agreement  for  both  types  of  samples  is  within  20%  of  theory  and  is  within  the  2-sigma 
error  bounds  of  the  data.  While  these  data  provide  the  necessary  calibration  and  confirmation  of  the 
theoretical  and  experimental  approach  taken  in  this  paper  further  refinements  of  the  measurements  and 
reduction  of  experimental  uncertainties  are  in  progress. 

In  Figure  7,  a  preliminary  residual  stress  map  for  a  polysilicon  cantilever  beam  is  presented  for  the 
critical  region  of  the  beam  near  the  support.  Locations  of  the  measurements  are  indicated  in  the  figure. 
Large  compressive  stresses  were  observed  despite  the  factor  of  2  uncertainty  that  arises  from  whether  the 
stress  is  uniaxial  or  biaxial.  The  measurements  were  taken  near  a  discontinuity  in  layered  films  of 
dissimilar  thermal  expansion  coefficients.  Further  study  is  required  to  resolve  the  stress  into  its 
components  for  this  geometry.  Preliminary  stress  profiles  in  terms  of  frequency  shifts  near  the  edge  of  a 
0.20-mm  square,  28-|xm  deep  trench  fabricated  using  the  Revise  Inc.  laser-enhanced  (thermal)  chlorine 
etch  are  shown  in  Figure  8.  The  square  trench  was  well  isolated,  greater  than  0.5  mm  from  other  trench 
structures  fabricated  on  the  same  wafer.  The  laser  has  a  spot  size  of  2-p,m  and  power  of  a  few  watts. 
While  some  variation  on  the  order  of  a  few  tenths  of  a  cm''  occurs  due  to  residual  stress  associated  with 
processing  defects  in  the  substrate,  the  fabrication  of  the  trench  has  caused  large  stresses.  Relative  to  the 
outside  of  the  trench,  a  large  negative  shift  of  1.9-cm''  is  observed  inside  the  trench.  This  indicates  that 
the  floor  of  the  trench  has  less  compressive  stress  (or  more  tensile). 

We  have  not  performed  the  measurements  to  determine  whether  the  shift  is  due  to  a  simple  uniaxial 
stress  or  some  combination  of  (biaxial)  stresses.  The  single  measurement  reported  is  also  insensitive  to 
shear  contributions.  The  simplest  assumption  is  the  case  of  only  uniaxial  stress.  More  work  is  needed  to 
fully  determine  the  contributions  from  all  of  the  full  stress  tensor  terms.  Assuming  this  is  due  only  to  a 
uniaxial  stress,  the  stress  difference  is  830-MPa.  The  magnitude  of  this  stress  is  sufficiently  large  to 
cause  significant  concern  for  the  mechanical  reliability  of  parts  using  this  fabricated  trench  approach.  In 
addition,  at  the  outer  and  inner  edges  of  the  trench  wall  sharp  stress  gradients  can  be  observed,  causing 
additional  concern.  The  origin  of  the  stresses  may  be  due  to  stress  concentration  caused  by  the 
discontinuity  of  the  edge  of  the  trench  or  to  heating  by  the  laser  to  cause  compressive  yielding.  Further 
experiments  are  required  to  identify  the  mechanism  responsible. 
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5.  SUMMARY 


In  summary,  we  have  demonstrated  analytically  that  the  two-dimensional  surface  stresses  normal  to  the 
<00 1>  plane  of  a  cubic  crystal  have  a  unique  relationship  with  the  Raman  frequency  shifts.  In  other 
words,  for  a  cubic  crystal  whose  spring  constants  p,  q,  r  are  known,  once  the  Raman  frequency  shifts  and 
polarization  characteristics  are  measured,  then  the  associated  surface  stress  components  can  be 
determined.  We  have  also  demonstrated  that  the  measured  Raman  frequency  shifts  for  a  single  crystal 
silicon  under  unidirectional  tensile  stress  in  both  (110)  and  (100)  directions  agree  well  with  analytical 
predictions.  With  this  validation  in  hand,  Raman  spectra  can  be  used  to  address  the  role  of  stress  for 
reliability  assessments  and  failure  mechanisms  in  MEMS  and  microelectronics.  The  use  of  these  Raman 
stress  diagnostics  will  lead  to  improved  lifetimes,  the  identification  of  defects,  and  improved  fabrication 
for  future  silicon  devices.  In  this  work,  high  stress  levels  were  observed  in  simple  single-  and  poly¬ 
crystalline  samples  that  raise  concern  for  deleterious  effects  for  these  types  of  structures. 

In  order  to  be  able  to  describe  the  stress  field  completely,  both  shear  stresses  and  normal  stresses  must  be 
measured.  At  this  time,  the  technique  for  measuring  shear  stresses  has  yet  to  be  developed,  although  the 
theoretical  background  for  doing  so  has  been  worked  out  (see  Reference  8).  The  Raman  measurements 
will  provide  the  maximum  benefits  in  stress  characterization  for  silicon-type  MEMS  systems  when  the 
entire  stress  field  can  be  derived  from  the  Raman  spectra. 


Figure  3.  Raman  spectra  for  unstressed,  single-crystal  Si  (data  collected  at  different  locations). 
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eq.  Shift,  R  Cl 


Figure  4.  Raman  spectra  for  polysilicon  deposited  on  Si02  on  single  crystal  Si  substrate  at  different 

locations. 


Figure  5.  Frequency  shift  as  a  function  of  external  uniaxial  applied  stress  in  the  (100)  direction. 
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Figure  6.  Frequency  shift  as  a  function  of  applied  external  uniaxial  stress  in  the  (110)  direction 


Figure  7.  Stress  map  of  a  20  x  200  |xm  polysilicon  (dashed  lines)  cantilever  beam. 


Figure  8.  Stress  profiles  near  the  edges  of  two  sample  specimens  of  a  0.20-mm  square  trench. 
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